Quantum Corrections to the Polarizability and Dephasing in Isolated Disordered Metals 
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We study the quantum corrections to the polarizability of isolated metallic mesoscopic systems using the 
loop-expansion in diffusive propagators. We show that the difference between connected (grand-canonical 
ensemble) and isolated (canonical ensemble) systems appears only in subleading terms of the expansion, and 
can be neglected if the frequency of the external field, <M, is of the order of (or even slightly smaller than) the 
mean level spacing, A. ff CO <C A, the two-loop correction becomes important. We calculate it by systematically 
evaluating the ballistic parts (the Hikami boxes) of the corresponding diagrams and exploiting electroneutrality. 
Our theory allows one to take into account a finite dephasing rate, y, generated by electron interactions, and 
it is complementary to the non-perturbative results obtained from a combination of Random Matrix Theory 
(RMT) and the o-model, valid at y — > 0. Remarkably, we find that the two-loop result for isolated systems with 
moderately weak dephasing, y ~ A, is similar to the result of the RMT-l-a-model even in the limit CO — !■ 0. For 
smaller y, we discuss the possibility to interpolate between the perturbative and the non-perturbative results. We 
compare our results for the temperature dependence of the polarizability of isolated rings to the experimental 
data of Deblock et al [Phys. Rev. Lett. 84, 5379 (2000); Phys. Rev. B 65, 075301 (2002)], and we argue that the 
elusive OD regime of dephasing might have manifested itself in the observed magneto-oscillations. Besides, we 
thoroughly discuss possible future measurements of the polarizability, which could aim to reveal the existence 
of OD dephasing and the role of the Pauli blocking at small temperatures. 



I. INTRODUCTION 

Interference phenomena in mesoscopic electronic systems 
require phase coherence, which is cut beyond the so-called de- 
phasing time T,j,. At low temperatures T < \K, where phonons 
are frozen out, dephasing is caused mainly by electron inter- 
actions, which lead to a finite dephasing ratd^J Y= 1/X(j). In 
large systems with a small Thouless energy, E-^h <C T, de- 
phasing crucially depends on dimensionality and geometry^ 
However, if the system is finite and T < Ej\^, spatial coordi- 
nates become unimportant and a OD regime of rather weak 
dephasing is expected to occui^. This regime is characterized 
by a universal temperature dependence of the dephasing rate, 
YoD '^ ^T'^ /E\^, where A is the mean-level spacing. This T- 
dependence of y can be explained by simple power counting: 
Pauli blocking restricts the number of available final scatter- 
ing states of the electrons, therefore both the energy transfer 
and the available phase-space are °^ T, similar to the standard 
result for a clean Fermi-liquid. However, despite of the funda- 
mental nature and the physical importance of OD dephasing, 
attempts to observe it experimentally in mesoscopic systems 
have been unsuccessful so far. 

In transport experiments, the OD regime is generally diffi- 
cult to observe, since quantum transport is almost insensitive 
to y at r ^ Ejh. For example, the weak localization correc- 
tion to the classical dc conductivity is cut mainly by the dwell 
time, Tdw ^ 1/Yod> see Ref. |4| for a detailed discussion. This 
is an unavoidable problem which occurs in any open system 
even if the coupling to leads is weak. 

In this work, we concentrate on interference phenomena in 
isolated systems, where Tdw -^ °° and where OD dephasing 
is not masked by the coupling to the environment. Deeply in 
the OD regime at y^ A, the spectrum of the isolated system is 



discrete^ and, in the absence of other mechanisms of dephas- 
ing, random matrix theory (RMT) can be used as a starting 
point for an effective low-energy theory at £ <C ^TtP^- Un- 
fortunately, RMT is not appropriate for a systematic account 
of dephasing. 

If one is interested in the (almost OD) regime y < A, where 
the spectrum is not yet discrete, the usual mesoscopic pertur- 
bation theory^ can be used, which is able to take into account 
dephasing in all regimes. However, the description of quan- 
tum effects in isolated systems provides a further technical 
challenge. Namely: the usual perturbation theory is well de- 
veloped for a fixed chemical potential /j, i.e. it describes sys- 
tems in the grand-canonical ensemble (GCE). Realizing the 
canonical ensemble (CE), where the number of particles is 
fixed instead, can be rather tricky, see, e.g., Ref. fl0|. In the 
following, we assume that a description in terms of the so- 
called Fermi-level pinning ensemble introduced in Ref. lITTI 
and 1 12 1 is applicable'^. 

The dephasing rate of an isolated mesoscopic system can 
be explored, for instance, by measuring quantum components 
of the electrical polarizability a at a given frequency (o: 



a(co)=d(co)-E(co)/|E((o)|^ 



(1) 



Here E is a spatially homogeneous electric field and d is the 
total induced dipole moment in the sample. 

Gorkov and Eliashberg studied the polarizability in the 
seminal work Ref. [14| by using results from RMT and 
found very large quantum corrections. Later, it was shown 
in Ref. ifTSl that the connections are significantly reduced if 
screening is taken into account correctly'^. Efetov reconsid- 
ered Gorkov and Eliashberg's calculation in Ref. 1 17 1 and de- 
rived a formula which accounts for screening in the random 
phase approximation (RPA) and expresses the quantum cor- 



rections to a in terms of correlation functions of the wave- 
functions and energy levels of the system. Subsequently, Mir- 
lin and Blantei** studied the polarizability using a combination 
of RMT and the diffusive o-model. In particular, they have 
calculated co-dependence of a at CO <C E-rh for the case of the 
CE at y = 0. Thus, neither the temperature nor the magnetic 
field dependence of a has been described until now. 

Besides the progress made in theory, experimental mea- 
surements of the quantum corrections have been reported in 
Ref [18l and [19|. The authors measured the T-dependence 
of the polarizability of small metallic rings placed in a super- 
conducting resonator (with a fixed frequency co) in a perpen- 
dicular magnetic field and tried to extract the T-dependence 
of T^ by using an empirical fitting equation. A fingerprint of 
OD dephasing was found at low temperatures, though a reli- 
able identification of the temperature dependence of T(p calls 
for a more rigorous theory. 

Motivated by the experimental results, we develop a pertur- 
bative theory for the quantum corrections Aa to the polariz- 
ability by using the mesoscopic "loop-expansion" in diffusons 
and Cooperons, where y plays the role of a Cooperon mass. 
We have chosen the experimentally relevant parameter range 
max(co,y) > A. Generically, the difference between the GCE 
and the CE can be important up to energies substantially ex- 
ceeding A, see the discussions in Ref. 1 10|. To check whether 
this statement also applies for Aa, we calculate leading and 
subleading corrections in the Fermi-level pinning ensemble. 
The former corresponds solely to the one-loop answer of the 
GCE while the latter includes the two-loop answer of the GCE 
and additional terms generated by fixing the number of parti- 
cles in the CE. We show that, within our approach, the lead- 
ing term of the perturbative expansion for Aa suffices for its 
theoretical description in the experimentally relevant param- 
eter range of Ref. ifTSl and |fT9l . This important result of the 
present paper allows us to find the dependence of Aa on tem- 
perature and on magnetic field. Our theoretical results are in 
good qualitative agreement with the experiments, though we 
show that the present experimental data is not sufficient for 
a reliable identification of OD dephasing. We suggest repeat- 
ing the experimental measurements with higher precision and 
lower frequencies and using the fitting procedures which we 
propose in the present paper. We have good hopes that the elu- 
sive OD regime of dephasing may be detectable in this manner 
in the near future. 

The rest of this paper is organized as follows: 

Section In] we derive a general expression for the polariz- 
abiliy as a functional of the density response function in the 
RPA. 



Section III we calculate the leading quantum corrections 
of the density response function for connected as well as iso- 
lated disordered metals. This part of the paper is rather for- 
mal and technical. Readers who are not interested in details 
of the calculations can safely skip it, paying attention only 
to our key results, which we list here. Firstly, we derive the 
one- and two-loop quantum corrections for the GCE which 
are presented in Eqs. ( 17|18 i of Subsection III A A "naive" 
loop-expansion for the GCE suffers from a double-counting 



problem of some diagrams which leads to a violation of the 
particle conservation law (electroneutrality) accompanied by 
artificial UV divergences. We suggest an algorithm of con- 
structing the diagrams which allows one to avoid all these 
problems. Our method can be straightforwardly checked for 
the one-loop calculations, see Fig. |2] and we extend it to the 
much more cumbersome two-loop diagrams shown in Fig. [3] 
Secondly, we calculate the leading diagrams which appear due 
to fixing the Fermi level in the CE. Their contribution is given 



by Eq. (|24| of Subsection pIB] 

Section ITVl we use the results from Section |lll] to derive a 
general equation for the quantum corrections Aa. 

Section IV] we compare our findings to the results obtained 
from a combination of the RMT and the a-model. We show 
that the diagrammatic result in the limit of a large conduc- 
tance, Eq. ( (30l l, qualitatively reproduces all features of the 
non-perturbative answers for almost OD systems at < CO < 
Exh, seeFig.|7] 

Section IVI] we apply our results for Aa to the ring geome- 
try, present a comparison with previous experiments and dis- 
cuss possible future measurements which can reliably confirm 
the existence of OD dephasing. 

II. POLARIZABILITY 

The polarizability ([T]) is governed by the induced dipole 
moment in the sample. 



d(co)= / d^x[x-nind(x,co)] 
Jv 



(2) 



where ni^d is the induced charge density. In the case of a good 
metal, screening should be taken into account in the random 
phase approximation (RPA), which results in the following 
expressions for the Fourier transform of rii^d'.^^ 



V ( fl frti 

e(q,co) 



(3) 



Here (|)ext(x,Co) = — E(co) x is the external electric potential, 
e(q,Co) ~ 1 — 2{/(q)x(q,K)) is the dielectric function, U is the 
bare Coulomb potential, and % is the density response function 
per spin. By using the Kubo formula, % can be expressed in 
terms of the commutator of the density operators n: 



x(q,«) 



d^x / df ([«(x,f),n(0,0)])e 



-i(qx-a)f) 



(4) 



We assume spatial homogeneity of the system, which is re- 
stored after disorder averaging. 

Inserting Eqs. ( 2]3 i in Eq. ([TJ, we find the following expres- 
sion for the polarizability: 



2e^ 



,x(q,«). 



«(to)= |^, .|2 y L 'l^ext(q,co) ^;;;'J (^ext(-q,co) . (5) 



q/0 



e(q,co) 



Note that the zero-mode does not contribute to a because of 
electroneutrality of the sample: 

X(qEEO,co)=0. (6) 



For a clean metal at co <?; vpq (vp is the Fermi velocity), % 
is local and is given by the density of states at the Fermi level: 



X(q,co^O)=Po. 



(7) 



The same equation holds true for a disordered (classical) metal 



at CO ^ Do^ {D is the diffusion constant), see Section [III 
Eqs.( |5|7| i yield the "classical" polarizability ao of the disor- 
dered sample. 



III. DENSITY RESPONSE FUNCTION 

In this section, we consider the density response function 
of the disordered metal which is needed to calculate the po- 
larizability, Eq. (|5]l. We will start with the loop-expansion of 
the disorder-averaged % in the GCE: xlAi=const = Z»- It is rel- 
evant for the polarizability of the connected system. Besides, 
the two-loop contribution to % is needed to study the differ- 
ence between the answers in the GCE and the CE. The latter 
is described in the second part of the present section. 

We consider only weakly-interacting disordered systems at 
small temperatures. The main role of the electron interaction 
is to generate a finite T-dependent dephasing rate for Cooper- 
ons. Therefore, we derive the density response function for 
the non-interacting system at T = and take into account 
y(r) at the end of the calculations. 



A. Grand canonical ensemble 

Simplifying Eq. Q for the non-interacting system at 
r = and fixed /j, x^ can be presented in terms of re- 
tarded/advanced {G^l^) Green's functions (GFs)^ 



X^(x,y,co) = -y de (8) 

p^+E(x,y)G^+g_„(y,x)+G^+E+„(x,y)p^+e(y,x) 



Here we have introduced the spectral function (or the non- 
local density of states): 



Pe(x,y) = ^[Gf(x,y)-G^(x,y) 



(9) 



In the presence of a random Gaussian white-noise disorder 
potential V (x) with correlation function 



y(x)y(y) 



1 



2jtpoT 

the disorder-averaged GFs are given by 

1 



§(x-y), 



Gf/^(k) 



e - Ek ± i/2x 



(10) 



(11) 



where T is the impurity scattering time and et is the particle 
dispersion relation. 




(b) 

^, {oy., {Aj^,, {sy.^ 



X .^ ■■ •■ 1^: 



FIG. 1. (a): One-loop correction to the density response function 
in the GCE. Retarded (advanced) GFs are denoted by solid (dashed) 
lines. Impurity lines, corresponding to the correlation function iflOl, 
are denoted by dotted crossed lines. Diffusive propagators are repre- 
sented by wavy double lines. They denote impurity ladders between 
the corresponding GFs of opposite retardation either in the particle- 
particle (Cooperon, Pc) or in the particle-hole (diffuson, P^) channel, 
(b) and (c): Dressed 4- and 6-point Hikami boxes which include di- 
agrams with one or two additional impurity lines connecting GFs of 
the the same retardation. 



The disorder average of Eq. (|8]l can be calculated with the 
help of the usual diagrammatic methods,-!- which yield the 
loop-expansion: 



X„(q,«)=Xo(q,co) 



-L5xiE- 



(12) 



Here j is the number of loops built from impurity ladder 
diagrams which include ladders in the particle-hole chan- 
nel (diffuson propagators) or in the particle-particle channel 
(Cooperon propagators). The leading (classical) term is well- 
knowrK 



Xo(q,co)==po 



Dq2 



Z)q2 - (CO 



(13) 



It obeys the fundamental requirement of electroneutrality, 
Eq. (|6]l, and reduces to Eq. (|7| at co <C Dc^ . 

The leading quantum correction SXgce describes the weak- 
localization correction to the diffusion constant^' and, there- 
fore, is also well-known. Nevertheless, we would like to re- 
call the basic steps of its derivation, which will be important 

to find the more complicated subleading term SXccp. 

The one-loop diagram, which yields §Xgce' i^ shown in 
Fig. flTa). It includes two diffuson propagators P^i and one 
Cooperon propagator P^, which are given by 



^rf(q,«) = 



Dq2 



na 



, fc(Q,co) = 



1 



DQ2-/(o + Y 



(14) 



The (ballistic) part of the diagram which connects the diffu- 
sive propagators is known as a 4-point Hikami box^^. It con- 
sists of three diagrams of the same order in (epx)^' shown 



in Fig. fltb) and labeled by {0}, {A}, and {B}, which are 
obtained by inserting additional impurity lines between GFs 
of the same retardation ("dressing" the Hikami box). The 
Hikami box should be calculated by expanding the GFs in 
each of the three diagrams in the transferred momenta and 
energies. A direct summation of the three diagrams gives 



>■ Q,£>) 



jj(direct sum) 
H4 



: 47tpoT'* [Dq2 +DQ2 - /co] . (15) 



The second and third terms in parentheses are manifestly in- 
correct as they violate electroneutrality, Eq. (|6|, and lead to 

an unphysical UV divergence in §Xgce- ^^^ incorrect terms 
originate from a double-counting problem: the diagram with 
a single impurity line, which contributes (via the diffuson) to 
the classical result of Eq. ([Tsj, is also included in the quan- 
tum correction ?>Xqc^ via the Cooperon attached to the "un- 
dressed" part of the Hikami box - the empty square {0}. One 
can eliminate unphysical UV divergent diagrams in the frame- 
work of the nonlinear a-model b y cho osing an appropriate 
parameterization of the matrix field22125j However, to the best 
of our knowledge a consistent procedure of their elimination 
in the framework of straightforward diagram techniques was 
not described in literature. As this is rather important for any 
calculation beyond the one-loop order, we give a detailed de- 
scription of such a procedure below. 

To avoid the double-counting, the Cooperon ladder of 
Fig.flTa) should start with two impurity lines when attached to 
the undressed box, while it should still start with one impurity 
line when attached to the dressed box. Thus, there is an am- 
biguity in the independent definition of the Hikami boxes and 
the ladder diagrams. We suggest a general algorithm which al- 
lows us to overcome this ambiguity and generate all properly 
dressed Hikami boxes obeying electroneutrality ^^ED. 

Let us consider the 4-point Hikami box shown in Fig. l2ta) 
to illustrate the method. Fig.l2ta) is obtained from Fig.flla) by 
"borrowing" two impurity lines to the undressed Hikami box 
from the attached Cooperon. We use this undressed box in 
Fig. l2ta) as a "skeleton diagram" which generates the dress- 
ings {A} and {B} of Fig. fllb) by moving one of the exter- 
nal vertices (with diffuson attached) past one of the borrowed 
impurity lines. Two possible movements of the left exter- 
nal vertex are indicated by arrows with lables {A} and {B} 
in Fig. l2|a). Fig. l2|b) shows all three components of the 
fully dressed Hikami box: two generated boxes, {A} and {B}, 
and the undressed box, {0}, where the external vertex is not 
moved. Dressing the Hikami box in this way removes the am- 
biguity, since all the Cooperon ladders attached to each of the 
boxes start with two impurity lines, thus avoiding the double- 
counting. Furthermore, using the identity^ 



Gf+a)(k+q)G^(k) 



q^O 



11 



1 - /xco 



Ge+co(k) - 



we illustrate in Fig. l2|c) that in the limit q — > the generated 
diagrams automatically cancel each other (to leading order in 
(Eft)^' ^ 1 ) at any Q and CO, thus ensuring electroneutrality 
and the absence of the UV divergence. 

Summing up the 3 diagrams drawn in Fig. l2lb) and us- 
ing the resulting expression to calculate the diagram shown 




{0} Q-k'/-..k' {A} Q-''/-..k' {B} k'+q/-.Q-k'-q 



"■■■X-"" ''■.Q-''-q k"+q X'--x- ■■■■' '\ 




-,Q-k-q 




^^^ Using Eq. (15): "^V "^^ -^ ( k 




smalt ' 



cancel cancel 

{0} + {A} + {B} = 0, for q ^ at any Q and (£>. 



FIG. 2. (a) The "skeleton diagram", which we use to generate the 
dressings {A} and {S} of the Hikami box shown in Fig.fTfb). The 
arrows with labels {A} and {B} indicate how the (diffuson attached) 
external vertex has to be moved to generate the correponding dressed 
boxes, (b) The resulting diagrams with the undressed, {0}, and two 
dressed boxed can be summed up directly, since no double-counting 
problem appears. To leading order in the transferred momenta and 
energies, (Dq^tjDQ^x, cot) <C 1, the sum of the three diagrams in (b) 
is 47tpT £)q . (c) Dressing the Hikami box by moving the external 
vertex guarantees that the answer vanishes at q — > 0, since the 3 dia- 
grams either cancel each other exactly (at any Q and to), or are small 
in this limit. This can be seen immediately after using the identity 
l |16^ and redrawing the boxes {0}, {A} and {B} as the 6 diagrams 
shown in the last line. 



in Fig. Ill a), we obtain the well-known resuli^ 



5XGCE(q'«) 






TZV (£)q2 - iCo)2 ■^ 



(17) 



Note that 5Xgce/Xo 



0(A/max(co,Y)), where A= l/(poV). 
Thus, Eq. ( [T7| i describes the dominating quantum correction 
to x^, if max(a),y) > A. 

To calculate the subleading quantum corrections, one has 
to consider the two-loop diagrams shown in Fig. [3] which 
contain momentum sums over diffuson or Cooperon propaga- 
tors, or both. Thus, their contribution is subleading in either 
G^(k)] (16) (^/inax(a),y)), (A/co), or (A/Dq^). Note that the diagrams 
J containing only diffusons are not relevant for the experiments, 

since they are magnetic field independent. We have used the 
algorithm described above to calculate the 4-point Hikami- 
boxes 7/4 ' of Fig. bI avoiding double-counting and main- 
taining electroneutrality, Eq. (|6]l. The "inner" Hikami box of 
^ is of different nature because it is connected 
Nevetheless, the same double 



Fig.[3j:g),//; 

to two internal Cooperons. 




(C) /'.(Qi.to) P,(Qj,(B) (d) /.,(Qj,a)) P,(Qi 




(2) 
FIG. 3. Diagrams contributing to the two-loop correction Xgce- ^^ 



swers for tiie Hikami boxes read: H^'^'^ 



r(dl,d2) 
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r(f2) 



0, w, 



H, 



m 



;ei) 

,(f4) 



xH. 



:<=2) 



(®(q(q + Qi + Q2)))^ Hi 

g2) 



X /f)'"' = 2©2(qQi)(qQ2), and /f|=^' = ©(Qj 



y/i*), see the main text for details; here 'D = 47ipx £). 



counting problem appears and can be overcome with the help 
of dressing this box by moving the vertices with the attached 
Cooperons. As a result, electroneutrality does not necessarily 

apply for i/4 ' , which is reflected by its y-dependence, see the 
next paragraph. Besides, the diagrams shown in Fig. [3tb-d) 
contain 6-point Hikami boxes. Their dressing is more subtle 
because of two issues, see the example shown in Fig.l4] which 
corresponds to the Hikami box H^ of Fig. pm): Firstly, start- 
ing with the undressed diagram and moving vertices into the 
attached diffusons, one cannot generate all required 15 dress- 
ings shown in Fig.[TJc). Instead, only 8 dressings can be ob- 
tained for the 6-point Hikami box, cf. Fig. |4| a). That problem 
can be solved by considering two more "skeleton diagrams" 
with one-. Fig. Wb), and two-. Fig. l4fc), additional impu- 
rity lines between GFs of the same retardation. All of the 
missing dressings can be obtained by applying the above de- 
scribed algorithm similar to Fig. Wa). Secondly, by moving 



(a) Q, Q, 

Q2 Qi 



{0'W^\ 



{-*'W'^\ 



r% *■■■/ '--S *■■■■ 



{B.oy 



:% 5t-/ 



{O.C}/ 



,,E \.> 



{-^•cv^'"s;\ {^-9"7^\ 



+ < 



:,^//^--:sJ/ 



■u // 



{"•WV^ ^-^^'^ 
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■<r\ ih-/ 



{B.D}rr\\ 



\\ +~; 




"fi-^ff'sX y^^ff'-^ 



'/...^ 5f,.> 






Qi Qi 



{CO}/ 7 



*■■ 



■ i 



+ 10 more diagrams 




V'0}£^'iA {-^Wi^^X (»-W% 
-X- / + '."C''^-.. .■ + ■■ ■^-. ) 



+ 14 more diagrams 



FIG. 4. Dressing of the 6-point Hikami box of Fig. [3lb) using the 
algorithm introduced in Fig. [2] (a) Only 8 of the 15 dressings are 
generated by moving the vertices. The other dressings are generated 
by adding one (b) or two (c) impurity lines, followed by repeating the 
procedure. This algorithm also generates products of 4-point Hikami 
boxes, indicated by a grey box. Summing up all 40 diagrams yields 
-12jrpT^£)q2. 



the vertices of the diagrams in Figs. |4tb,c) new diagrams of 
the same order in (epT)^' ^ 1 are generated, which look like 
products of two dressed or undressed 4-point Hikami boxes 
with a few-impurity ladder in-between. Several examples are 
highlighted by grey boxes in Figs. l4tb,c). It is not a priori 
clear whether such diagrams belong to the diagram shown in 
Fig.[3lb) or Fig.lSje). However, keeping them only in the di- 
agram Fig. |3jb) allows us to maintain the electroneutrality in 
all two-loop diagrams. The total result for H^ is obtained by 
summing 40 generated diagrams. The 6-point Hikami boxes 
of Figs.[3lc,d) can be calculated analogously. 

Before presenting the final answer, we would like to dis- 
cuss the way of how to reinstate the finite dephasing rate in 
the equations. Firstly, y must be included as a mass term 
in all Cooperon propagators. Secondly, when calculating the 

Hikami box H^'" of Fig. Pig), only the number of coherent 
modes has to be conserved. The latter is in contrast to all 
other Hikami boxes, which obey the usual electroneutrality 
condition, i.e., the conservation of the total number of parti- 

(g2) 

cles. Hence, H^ is the only Hikami box of the two-loop 
calculations which is sensitive to dephasing of the Cooperons. 
This statement can be checked directly with the help of the 



model of magnetic impurities. Introducing a slightly reduced 
scattering rate for all elastic collisions in the particle-particle 
channel, 1/t— > 1/t — Ymi, where Y^i ^ l/x, and keeping l/x 
for collisions in the particle-hole channel, we observe that 
the Cooperon acquires the mass ymi since magnetic scattering 
breaks time-reversal symmetry. Hence, magnetic scattering 
rate is similar to the dephasing rate; they both provide consis- 
tent infrared cut-offs for Cooperons. Applying the algorithm 
described above, we find that, among all the two-loop dia- 
grams in Fig. [3] the rate Ymi appears only in the expressions 
for Cooperons and in the Hikami box H^ . In the latter case, 
it leads to changing DQi^ to DQi^ + Ymi- Using the analogy 
between magnetic scattering and dephasing, we conclude that 

Y enters H^ in the same way. 

Omitting lengthy and tedious algebra which will be pub- 
lished elsewhere, together with a detailed proof of the validity 
of our method and an analysis of the IR cut-off in systems 
with magnetic impurities, the answer for SXq^e reads: 

1 2mDq^ ^ 



5Xgc£ (<!'«) 



P,(Qi,co)P,(Q2,co) 



(27l)2poy2 (Dq2 
/Dq2 



iCO 



4D(qQi)(qQ2)/q^ 

/CO 



(18) 



V£>q2-/co Z)(q + Qi + Q2)2 

MQi,mi{Q2,^) [ z)(q + Qi+Q2)2-/co 
-i'c(Qi,«)P,,(Q2,co) (2 + 2/coP,(Qi,co)) 



-1 



To conclude this section, we would like to note that our 
method of dressing the Hikami boxes goes far beyond the 
initial ideas of Ref. |24|. It is a very powerful and generic 
working tool which can be extended to even more compli- 
cated diagrams, including higher loop corrections, and non- 
trivial physical problems. For example, our method can be 
straightforwardly used to describe mesoscopic systems in the 
ballistic regime, cf. Ref. 1261 . Therefore, the diagrammatic 
approach presented above is complimentary to the diffusive 
nonlinear a-model which fails to yield ballistic results. One 
can invent alternative digrammatic tricks which help to avoid 
the complexity of the Hikami boxes with scalar vertices. For 
instance, the density response function can be obtained by cal- 
culating the current response function (averaged conductivity) 
first and then using the continuity equation. In the latter ap- 
proach, the dressed scalar vertices are replaced by undressed 
vector ones, which greatly simplifies the calculation^^. How- 
ever, this method cannot describe the full q-dependence of %, 
which is crucial for the polarizability. We have checked that 
both approaches give the same results in the small-q limit. 



B. Canonical ensemble 

In this section, we study the disorder average of the density 
response function % in the CE, where the number of particles 
A^ is fixed in each sample. Let us first discuss the properties of 
the statistical ensemble which corresponds to the experimen- 
tal measurements of the polarizability, such as the experiment 



discussed in Section VI We are mainly interested in the be- 
havior close to the OD regime, where due to x^, > l/Ej^, there 
is no self-averaging. Instead, the disorder average is usually 
realized by an ensemble average. The samples from the en- 
semble differ in impurity configuration and can have slightly 
different particle number At T = (in the ground state) all 
single-particle levels below the Fermi level 8f are occupied. 
However, one cannot fix Ef for the whole ensemble due to 
randomness of the energy levels and due to the fluctuations 
of A^ from sample to sample. This can be taken into account 
by introducing an Ep which fluctuates around the typical value 
/j"," /j*^ fixes the mean value of A^ in the entire ensemble. 
It has been shown that such ensembles of isolated disordered 
samples with fluctuating Ef ca n be d escribed by the so-called 
Fermi -level pinning ensemble, '^ '^ which is realized as fol- 
lows: (i) the Fermi-energy is pinned to an energy level Ek, 
such that Ef = Ek + 0. (ii) the level Ek is sampled from a 
weight function P(Ek), which is centered at /j" and is nor- 
malized: /P(E)dE = 1. The support of P(Ek) should be much 
smaller than /j" but much larger than A. The correlations re- 
sulting from fixing A^ in the given sample are subsequently re- 
duced to the additional correlations induced by disorder with 
the help of the following procedure: The expression for the 
density response function averaged over the fluctuating Fermi 
energies and over disorder reads: 



1 



Lk^i^k) k 



(19) 



In Eq. ( [T9] l we have assumed that the numerator and denom- 
inator can be averaged over disorder independently, see the 
discussion in Ref. 1111 . Since the averaged density of states 
depends only weakly on disordei^l and is almost constant on 
the support of P, the denominator of Eq. ([T9| can be simpli- 
fied 



^P(Ek)=y/ dEP{E)p^^poV. 



(20) 



Inserting Eq. ^ and Eq. ( pOJ i into Eq. ([19]), we find the disor- 
der averaged density response function in the CE: 



5c(q,co) = -/ d£F(£)p£Xi-(q,co)-5C^ + 85Cc£- (21) 
Po J -00 

The loop-expansion of % was calculated in the previous sec- 
tion. The quantity 5xc£ describes additional contributions 
resulting from fluctuations of Ef. It is governed by the irre- 
ducible part of the integrand: 



Po 

In Eq. ( |22| i, we have assumed that the disorder-averaged quan- 
tities are (almost) independent of the absolute values of the 
particle energies. As a result, the exact form of the weight 
function P(Ek) is not important. Let us now derive the leading 
contribution to 5%c£- 



§Xc£ = — / dEP{E) (p£X£(q,co)-p£X£(q,co)^ 

— (pirX£(q,«)-poX£(q,«)) . (22) 
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FIG. 5. Diagrammatic representation of the term pEXEi^.,^) from 
Eq. \12) before impurity averaging, cf. Eqs. l|8|9[). 



Eq. ( [T8[ l. As a result, yhas to be included only as a mass term 
in the connected Cooperon. 

Summing up all parts and calculating the auxiUary deriva- 
tive, Eq. (|23]l, we obtain the one-loop answer for ?>Xce'- 



(2^)^PoV^L? ^'"^ 



5xc£(q,«) 



Oq-^-ico 



^a(Q,co)P„(Q,0)+Pa(Q + q,ra)Fcx(Q,0) 



(7c/j(Q,-e+ra+X) 
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FIG. 6. One-loop diagrams whicli contribute to tfie disorder averaged 
5xc£i Eq. l[22}, before taking the derivative 3/3A,, cf. Eq. ( |23| l. Both 
4-point Hikami boxes in (a) are given by /f^ = 27lpot (£)q — to) . 



Electroneutrality is restored in Eq. ( p4| i after summing all the 
diagrams of Fig. [6] Thus, all contributions, Eqs. ( 17 1, (flSll and 
(24 1, obey the electroneutrality condition; therefore, %(q = 
0,C0)=0. 

Note that the one-loop contribution Sx^i^ ( |24| i is of the 
same order in (A/max(a),Y)), (A/co) or (A/Dq^) as the two- 

(2) 



Eq. (flSll. That is why ht^L is 



loop contribution Sx^^^, 

needed to analyze the difference between the GCE and the 

CE for energies of the order of 0(A). In the following, we 

will often refer to ?>Xqce as the result from "1st order" pertur- 
bation theory, and 5Xgc£ + ^gce + ^^ce (or ^^gce + ^gce') 
as the result from "2nd order' perturbation theory for isolated 
(or connected) systems. 



Diagrammatically, the additional factor p£ in Eq. ( [22| l is 
represented as a closed fermionic loop with a vertex between 
two (disorder averaged in further calculations) GFs which 
have the same retardation, energy and momentum, see Fig. E\ 
Following Ref. 1.28 J , we greatly reduce the number of possible 
diagrams in Eq. ( |22) l by generating this vertex with the help of 
an additional energy derivative: 



Gf^(k)2 



,. —R/A 



(k). 



(23) 



IV. QUANTUM CORRECTIONS TO THE 
POLARIZABILITY 



The quantum corrections to a can be found after inserting 
the decomposition X = Xo + 8x into Eq. (|5]l and expanding 
the density response function in the RPA, x/e, in 8^. Note 

that the latter can contain Sx^j^i and Sx^^ depending on 
the ensemble which we consider and on the accuracy of the 
loop-expansion. This expansion up to terms of order (9(5x)^ 
yields: 



After disorder averaging, we find two types of one-loop 
diagrams which contribute to 5xce, see Fig. [6] (i) the di- 
agrams in Fig. |6|a) are obtained by pairing the closed loop 
with the G^G^ terms of Xe (first term of the second line of 
Fig. B); (ii) the diagrams of Fig. l6tb) result form pairing with 
the G^G^IG^G^ terms (second and third term). Fuithermore, 
4 more diagrams can be constructed where Cooperon propa- 
gators are replaced by diffuson ones. 

The double-counting problem does not appear in the di- 
agrams in Fig. l6la), which contain 4-point Hikami boxes. 
Therefore, the method which we used for the GCE diagrams 
is not needed here. The only subtle issue in their calculation is 
that the diagrams are small if the closed loop, p£, is connected 
to the bubble, Xfi. by only one single impurity line. Thus, at 
least two such connections must be taken into account either 
in the ladder (which starts then from two impurities) or in the 
ladder (which can start from one impuiity) and the particular 
dressing of the Hikami box which connects p£ to X£. Further- 
more, the 4-point Hikami box in Fig. [6|a) does not acquire a 
dependence on dephasing rate y, which can be checked with 
the help of the model of magnetic impurities discussed before 



x(q,to) _ xo(q,to) 

e(q,co) "^ eo(q,co) 



1 



8x(q,«) 



eo(q,co)xo(q,co) 

l-£Q(q,(o) / 5x(q,«) 
eo(q,(o)2 VXo(q,(o) 



(25) 



where eo(q,co) = 1 — 2C/(q)xo(q,K))- To separate the fre- 
quency dependence due to classical diffusive screening from 
the frequency dependence of the quantum corrections, it is 
convenient to rewrite Eq. (|25| as follows: 



x(q,«) 

e(q,co) 



po5(q,co) 



.-*^n„». 



(26) 



+ 8f/(q)xo(q,co)4^^F(q,co)2 

^(|q| ) 



Here we have introduced two dimensionless functions: 

Dq^ 



5(q,(o)=l-2f/(q)po- 



(27) 



which describes classical diffusive screening, and 



— 2nd order pert, theory (y = A/7i) — RMT+a-model 



F(q,(0) 



(£)q2 - iaf 
Z)q2 



7ty5x(q,co), 



(28) 



which describes the quantum corrections to %. g(L) denotes 
the dimensionless conductance of a diffusive system of size L: 



g{L) = 27i£Th(i)/A , E-niL) = D/L^ 



(29) 



Eqs. (|26|-(|28]l together with Eqs. ([T7]i, ([18]) and (|24]) are the 
first major results of this paper. The quantum corrections 
Aa are obtained by substituting the terms ~ F and ^ F^ of 
Eq. ( p6l l into Eq. (|5]) and summing over q. We remind the 
reader that the zero mode does not contribute to the polariz- 
ability due to electroneutrality x{<d,(i)) = and, therefore, we 
can assume |q| 7^ in Eq. p6l ). The typical momenta which 
govern the sum in Eq. (|5]l are |q| ^ \/L since the external po- 
tential (|)ext varies on the scale of the sample size L. But we 
will keep q below for generality. 



V. COMPARISON TO RMT+o-MODEL 
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Let us now compare the results of our perturbative calcula- 
tions with those of Ref.|8| which are obtained from a combi- 
nation of the RMT approach and the nonlinear a-model. The 
latter will be referred to as "RMT-na-model". This compari- 
son requires an assumption Ej\^{L) ^ max(A,a),Y) which in 
particular means g{L) — > °°. In this limit, the term ^ F^ in 
Eq. ( |26l l acquires an additional smallness (which can be es- 
timated as 0{\ / g{L))) and can be neglected while the term 
~ F^ becomes independent of q. Next, we keep only the 
zero mode contributions in all sums over internal momenta 
in the expressions for %q'(jI and Sx^g and consider the dif- 
ference of F calculated for unitary and orthogonal ensembles: 
dBF{&) — F{(i),B — > 00) — i^(co^O), where B is the sti-ength of 
an external magnetic field. The terms which contain only dif- 
fusons are cancelled in ?>bF. 

Using Eqs. ([T7]i, ([T8]l and (|24]), we obtain 



/CO A 



(30) 



y-Zco 2k 

"X-GCE 



2y 



/CO - 2y 

(y-/co)2 7(7- /co) 



"f-acE 



^l 



Subscripts under the braces explain the origin of the corre- 
sponding terms. The last term must be taken into account only 
in the CE. The counterpart of Eq. ([30]l obtained from RMTh-o- 
model in Ref.|8 1 reads: 



RMT + a: 5bF(co) = 1 



Me / 1 



1 



A Ve-K> e + co 



CE 



/■E-O 

e5B;?2(e)+A5B/?2(e) + y dei 85^3 (£,ei) 



GCE 



FIG. 7. (color online) The quantum corrections to the polarizability 
in the limit £'Th(^) 2> max(A,co,Y) for the GCE (upper panel) and 
the CE (lower panel). We compare real (solid lines) and imaginary 
(dashed lines) parts the function 8bF obtained from 2nd order per- 
turbation theory, Eq. l |30| l, and from the RMT+o-model, Eq. ( |31[ l. 



Here 7?2,3 are the usual (dimensionless) two- and three-level 
spectral coiTelation functions, ^3(e,ei) — 7?3(e,ei) — 7?2(e), 
and 8b7?2,3 denotes the difference of the coiTelation functions 
without and with time-reversal symmetry. We have marked in 
Eq. ( |3T| ) the relevance of different terms for the GCE and the 
CE. 

We remind the reader that the RMT-na-model results are 
valid for y = and cannot straightforwardly describe a y- 
dependence, while our perturbative result, Eq. ([30|, is valid 
only if A < max (y, co) . To resolve this issue, one should set in 
Eq. (|30]l y ~ A. Eq. (|30]) yields 5bF(co -^ Q,g -^ °°) = A/(7uy) 
for the GCE. Therefore, we have chosen y = A/k to ensure 
the correct limit 5bF(co — > 0,g — > °°) Ip^E ^ ^■ 

The comparison of the results obtained from RMTh-o- 
model and from the perturbative calculations are shown in 
Fig. [7] for the GCE and the CE. Apart from the oscillations 
in the RMT-no-curves, whose origin is non-perturbative, the 
agreement is excellent. The asymptotic limits are fully re- 
covered in the perturbative calculations: (i) 5bF{(x) ^ A,g —> 
00) — > 1 for the both ensembles; (ii) 5bF(co — ?► 0, g — > 0°) — > Q 



(2) 



(1) 



/oi\ in the CE due to cancellation of Sx^^^ and Sx^^- The latter 
property of the CE holds true at any y in 1st and 2nd order 
peiturbation theory. 

We conclude this section by noting that the perturbation 
theory is able to reproduce the results of the RMT-no-model 
with good qualitative agreement, which is the second major 
result of our work. 



VI. POLARIZABILITY OF AN ENSEMBLE OF RINGS 

The experiments described in Ref. ifTSJI and lfT9l were done 
on a large number of disordered metallic rings. The rings 
were etched on a 2D substrate and were placed on the capac- 
itative part of a superconducting resonator, where a spatially 
homogeneous in-plane electric field E(a)) acted on them. In 
terms of the coordinate along the ring, x G [0,2jti?], where R 
is the ring radius, the external electric potential of this field 

is (|)ext(ji;,H)) — |E(a))|7?cos(jc//?) +(|)ext' ^^d its Fourier trans- 
form reads 



(j)ext(?,CO) = -|E(C0)|/?27I [?>^AlR + ^q,-llR] 



-€. 



Jq.O ■ 



k(0) 



(32) 



The constant shift of the potential (|)gxt does not contribute to 
the polarizability. Therefore, the sum in Eq. (|5j involves only 
two modes, q = 1/R and q = —l/R, which yield 



a(a)) 



4e2 



|E(co)p27i/? 



^ext( 



,co) 



X(?,«) 



£(?,«) 



-l/R 



= 2.e^-R^^^"''^^^ 



(33) 



e(l/i?,co) ■ 

In Eq. ([33|, we have taken into account the symmetry of the 
summand under the inversion q — > —q. 

The Coulomb potential in quasi- ID is given by 

U{q)^2e^\n{\qW\), |?W|<1; (34) 



where W <^R is the width of the ring. Inserting Eq. ([34J into 
Eq. ( |27] i, we find the screening function of the quasi- ID ring 

at ^= l/R: 



S{l/R,K>)= l + {KW)\n{R/W)/% 



1(0 



-) 



-1 



(35) 



kW>1 



{KW)ln{R/W) 



5o < 1 . (36) 



We have introduced the 2D Thomas-Fermi screening vector, 
K — 4ne^po/W with po being the quasi-lD density of states, 
see e.g. Ref. f9l, and assumed sufficiently strong screening, 
kW ^ 1, such that S reduces to the co-independent constant 
^o. This agrees with the experiment where one can estimate 
(KW)ln(i?/W) « 18. Therefore, we focus below only on the 
limit of strong screening. Note that in this limit, the product 
U {I / R,(}li)S[\ / R,(£i) can be also simplified 



f/(l//?,co)5(l//?,co)«-l/2po. 



(37) 



The classical part of the polarizability comes from inserting 
the leading term of the expansion (|26]l into Eq. ([33]): 



ao ~ iKe^R^ PqSq 



nR' 



2\n{R/W) 



(38) 



Using Eqs. ( 13|37 i in Eq. ( p6l ), and inserting the result into 
Eq. ( |33| l, we obtain the quantum corrections to the polariz- 
ability up to the term ^ [F /gf-; 



Aa(co) F{R-\(£,) ^ ETh{R) fF{R-\oi) 



25oao 



8iR) 



-2 



£ThW-/coV 8{R) 



(39) 



0) = 0.lA,g = InEn/A = 10 
1st order pert, theory: 2nd order pert, theory: 

^^ Re5gf''' Re5fif<''+Re5BFP) 

■■■■■ g^°°[Eq.(30)] g^°° [Eq.(30)] 




FIG. 8. (color online) Comparison of perturbative 1st order, 5bF''', 
2nd order, 5bF(') +5bFP', and interpolated (to the RMT+a-model 
limit) results for the quantum corrections to the polarizability in the 
parameter range CO < A < Ej\^. 



Let us regroup the terms in Eq. ( |39| l to single out the terms of 
1st and 2nd order perturbation theory: 



Aa(co) ^ 1 



-(f(i)(1/7;,co) + fP)(i//;,co)) (40) 



with 



F(^)(l//^,co) = (2.^/^)i^Ilg^5x[;^,(l//^,co), (41) 



and 



F^^\l/R,&) 



(42) 



(2.-/^) (^^g)^^^'")' (5x[?^,(l/^,co) + 5x[:^(l/^,co) 
-^(2.^/.)^(^lfci^(5xg^,(l//.,co)"^ 



We emphasize that all three parts of the density response 

function, ?>%qq^ and Sx)-.^, are generically important for the 
theoretical description of the experimental data with the help 
of Eq. ( |40]i if the rings are isolated. Having obtained Eqs. ^V7\ , 
( 18]) and"p4| (and Eq. ([30]) for the limit g^°°) and Eqs. ( p9p - 
( 42 1, we are now in the position to analyze different options to 
fit the experimental data. Ref. ifTSl and |fT9l focused on the T- 
dependence of the real part of the quantum corrections, thus, 
in the following we will concentrate on ReAa. 

The crossover to OD dephasing occurs when y decreases be- 
low A. We expect that the ideal parameter range to study this 
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solid: ReS^Ffi), dashed: ReS^ (fC'+F^^') , dotted: j^J^j^ 




y=0.1A 
T=0.5A — 
T=1.0A — 



FIG. 9. (color online) Comparison of perturbative 1st order, 8bF K 
2nd order, 5bF('' +8^/^^', and interpolated (to the RMT+o-model 
limit) results for the quantum corrections to the polarizability (a) as 
a function of y for different values of CO and (b) as a function of CO for 
different values of y- 



crossover experimentally in the CE is (O < A < Eji^. How- 
ever, it is important that the conductance should be only mod- 
erately large, since Aa is suppressed in the case of extremely 
large g, cf. Eq. ( |39| l; and the frequency should not be too 
small, since the quantum coiTections to the polarizability of 
isolated systems are suppressed in the static limit, see Fig. IT] 
Let us first discuss our general expectations for this parameter 
range, which are illustrated in Fig. [8] The simplest regime 
is 1 ^ Y/^ ^ 8 where the loop-expansion can be justified 
and the difference between the GCE and the CE is negli- 
gible. Keeping only the leading term, we obtain a power 
law for the dependence of Aa on y. This power law can 
be derived straightforwardly after noting that, in the range 
(y, a))/A ^ g, one can use the approximation Eq. ( pO] ) and 
find ReAa ^ Re5BF(co) ^ co^/y^ for co <C y. 

The subleading terms, which in particular describe the dif- 
ference of the GCE and the CE, are able to improve the theo- 
retical answer for y being slightly smaller than A. However, 
ReSflF^^' (and, correspondingly, the difference between the 
ensembles) is small at any y for moderately small frequen- 
cies, see the example co = 0.4A in Fig.M Therefore, 5bF^^' 
suffices to fit the experiment at CO > 0.4A. The T-dependence 
of Aa saturates to the value predicted by the RMT+o-model 
at Y < CO which makes the range of pronounced OD dephas- 
ing (co < Y ^ A) too narrow even at CO ~ 0.4A, thus, smaller 
frequencies are needed. Of course, the perturbation theory 



is no longer valid if both co and y are small. In particular, 
when F'^' becomes of order of F''' it can lead to changing 
the overall sign of RebsiF^^^ +F'^^'), see the cut of the lines 
in Fig.l8]marked "pert, theory breaks down". We believe that 
this sign change is unphysical and, moreover, it contradicts 
the prediction of the RMT-na-model. Nevertheless, our calcu- 
lations show that the power law, which is obtained in the per- 
turbative region from the leading correction, can be extended 
well into the non-perturbative region CO < y ^ A. This pro- 
vides us with the unique possibility to detect the crossover to 
OD dephasing directly from the amplitude of Aa. It is in sharp 
contrast to the quantum corre ctions to the conductivity, which 
always saturate at y < A.!^!^ 

Let us illustrate our unexpected statement with the help of 
Fig. Is] We know the exact value of Aa in the limit y ^ 
from the RMT-no-model and the correct behavior of Aa for 
y being of order of (and slightly below) A. Using these ref- 
erence points, one can interpolate the dependence 5a(y) for 
the whole region < y < A. Since the slope of the interpo- 
lated curve is only slightly different from the perturbative one 
for y > 0.3A, the leading answer of perturbation theory can 
be used to detect the crossover to OD dephasing. If the range 
y > 0.3A is not sufficient for unambiguously fitting the exper- 
iment, the whole interpolated curve can be used instead. 

The authors of Ref. ifTSll and lfT9ll used a superconducting 
resonator with fixed frequency CO ~ 0.2A ~ 17mK to measure 
Aa(r) of the rings. In the following we will apply our theory 
to explain the experimental results of these papers. We note 
that the qualitative difference in the slope of the curves ob- 
tained from the three options for fitting - (i) the interpolated 
curve, (ii) the result of 2nd order perturbation theory, and (iii) 
the leading perturbative result - becomes rather insignificant 
at CO ~ 0.2A and y > 0.3A, see Fig.l9ta). The main difference 
between (i) and (iii) is that the saturation originates at slightly 
larger y than the leading perturbative result would suggest. 
Thus we can safely keep F^^' and neglect F'^' to fit the data, 
which makes our task simpler^^. 

The experimental results for the ring polarizability can be 
distorted because of a parasitic contribution from the res- 
onator The latter has been filtered out in the experiment with 
the help of an additional weak magnetic field B applied per- 
pendicular to the rings, such that Aa becomes a periodic func- 
tion of the magnetic flux through the ling. Measuring the T- 
dependence of the (|)o/2 oscillations, cf. Fig. 9 of Ref. |19|, 
allows one to focus purely on the response of the rings. Using 
Eq. ([17]) in Eq. (|28]), we find 



Aaoc f(^'(1/7;,co) = /co^P,.(Q,(o) 
Q 



(43) 



= /COL /°°df e'™ y _^e-(«^)V40rg.e«^-Y, 
Jo V V47t£)? 

where 9 = 47t(|)/(|)o and (|) is the flux through one ring, and 
L = 2kR. Taking the Fourier transform and selecting the (|)o/2 
signal gives: 



5w2^^'^(l/^,«) = 



iCO exp (^- y/{y-i(a)/ETh{L) 
^Eji,{L){j-i(o) 



(44) 
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(a) Re5A„/2F(') CO = 0.2A g = 27t£Th/A = 5.6 
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(b) Measured amplitude of the isfo/l oscillations (Ref.[20] and [21]) 
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FIG. 10. Amplitude of the (])()/2 oscillations, (a) Expected de- 
pendence on Y from our theory, Eq. \A4\ , for the parameter range 
CO <C A <§; iijh. (b) Experimentally measured data as a function of 
temperature and possible interpretation. Note that the theory (see e.g. 
Ref. l33l and (34)) predicts Yod °= T^ in the OD regime, and Yerg °= T 
in the ergodic regime, therefore, the y^ ' behavior indicated in (a) 
encompasses both the T^^ and T^^l'^ behavior seen in (b). 



The function 5, 
5bF(i), cf. 



in o,j,jj/2/^'''' is shown in Fig. 10 a). It is similar to 
Fig. 8 however, the dependence of h^^nF^"^' 



on y 



is governed by a <=<: f^l'^ power law in the regime co <C Y<C ^A, 
and in the regime gA <C Y, the (|)o/2 oscillations are exponen- 
tially suppressed. The theory predicts a OD dephasing rate, 
Yod = flAr^/Zilj^, at low temperatures and an ergodic dephas- 
ing rate, Yerg = bAT /Ej^t, at higher temperatures, where a and 
b are system-specific, dimensionless coefficients of order ^ 1, 
see Ref. 1 33 1 and 1 34| . The crossover between the two regimes 
occurs at a temperature Tiross = a^^Th- We expect that the sat- 
uration at Y = CO occurs in the OD regime, corresponding to a 
temperature T^at = -^-£'Th-\/K'/A. Note that the conductance 

of each ring was rather small, g{L) w 5.6, such that the Thou- 
less energy EjilL) « 0.9A. Thus, depending on the coeffi- 
cients a and b, Tdoss and Tsat can be relatively close to each 
other. 

The experimental result for the T dependence of the (|)o/2 
oscillations is shown in Fig. flOfb). The measurements were 
done in the temperature interval CO ~ 0.2A <T< 4A. Based 
on the preceding discussion, we offer the following interpre- 
tation of the data: At low temperatures T < 1 .2A, the quan- 
tum corrections depend only weakly on T and are almost sat- 



urated. At intermediate temperatures 1.2A < T < 2.5A the 
slope of the data is steep and consistent with OD dephasing 

Aa(r) °c Yod''^ °^ ^ "*■ ^^ higher temperatures T > 2.5A, the 
slope of Aa(r) decreases and is consistent with ergodic de- 
phasing Aa o<: Yerg °= T^^'^. The crossover temperatures, 
Tsat — 1.2A and Tcross — 2.5A, correspond to coefficients a ~ 
0.1 and b ~ 0.3, which are close to the values predicted in 
Ref 1 33 1 (a ~ 0.04 and fe ~ 1). However, we stress that this 
interpretation is based only on very few data points, and we 
do not claim that the experiment clearly shows a crossover to 
OD dephasing. Further experiments are needed to support this 
statement, see Section I VIII 



VII. CONCLUSIONS 

Understanding interference phenomena and dephasing in 
mesoscopic systems at very low temperatures is a subtle 
issue which has provoked controversies between different 
theoretical approacheg^^, as well as between theory and 
experiments^^. Quantum transport experiments cannot give 
a certain answer to all questions because of unavoidable dis- 
tortions due to the coupling to the environment. The response 
of isolated disordered samples, on the other hand, provides a 
"cleaner" setup to study dephasing, and gives one the possi- 
bility to settle long-lasting open questions. 

We have studied the quantum corrections to the polariz- 
ability of isolated disordered metallic samples aiming to im- 
prove the explanation of previous experiments (Ref. ITSl and 
fT9\), and to suggest new measurements, where the elusive 
OD regime of dephasing can be ultimately detected. Using the 
standard strategy of mesoscopic perturbation theory, i.e. the 
loop-expansion in diffusons and Cooperons, we have devel- 
oped a theory, which (i) accounts for the difference between 
connected (GCE) and isolated (CE) systems, and (ii) is able 
to describe the low frequency response of disordered metals, 
taking into consideration weak dephasing induced by electron 
interactions. We have shown that the difference between the 
GCE and the CE appears only in the subleading terms, there- 
fore, we have extended the calculations up to the second loop. 
An important by-product of these calculations is a system- 
atic procedure to evaluate the Hikami boxes, see Fig. l2] and 
|4J which is based on a fundamental conservation law^'^: elec- 
troneutrality of the density response function. Our main ana- 
lytical results for the quantum corrections to the polarizability 



are presented in Eqs. (26i-(28 i with Eqs. (17i, ( 18 1 and (24 1. 



We have demonstrated that, in the experimentally relevant 
parameter range, the difference between the statistical ensem- 
bles is unimportant and one can fit the measurements by using 
the leading term of the perturbation theory. The authors of 
Ref. HSl and Ull have tried to find OD dephasing with the 
help of an empirical fitting formula. By using the more rig- 
orous and reliable Eq. ( |44] i, we have confirmed that OD de- 
phasing might have manifested itself in the T-dependence of 
magneto-oscillations at T < Ej^. Unfortunately, the T-range 
of interest here is rather narrow, and only few experimental 
data points are available there. Therefore, we are unable to 
claim conclusively that OD dephasing has been observed in 
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the experiments. However, we can straightforwardly suggest 
several experiments which might yield conclusive evidence 
of OD dephasing: Firstly, one can repeat the measurement of 
Ref 1 18 1 and |19|, but with a larger number of data points 
around the crossover temperature Tcross, see Fig. 10 while si- 
multaneously improving the measurements precision. Since 
the theory predicts a drastic increase in slope of the (|)o/2- 
oscillations at the crossover (from T^^'^ to T^^), even such 
measurements should be able to reliably confirm the existence 
of OD dephasing, thereby uncovering the role of the Pauli 
blocking at low T . Secondly, it is highly desirable to extend 
the T-range where the crossover to OD dephasing is expected 
to appear, which can be achieved by decreasing co and/or in- 
creasing g. However, a very large conductance and ultra-small 
frequencies are nevertheless undesirable, because in these lim- 
its the quantum corrections to the polarizability are reduced. 
Thus, improving the precision of the measurement is needed 
anyway. Besides, fitting with the help of the leading pertur- 
bative result fails at very small frequencies, see Fig. l9] This 



difficulty can be overcome by taking into account our two- 
loop results and/or using an interpolation to the y — > limit 
from the RMT-no-model, see Fig. [8] 

To summarize, we have shown that the quantum corrections 
to the polarizability are an ideal candidate to study dephasing 
at low T and the crossover to OD dephasing. We very much 
hope that our theoretical results will stimulate new measure- 
ments in this direction. 
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